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Abstract 



The spectrum and eigenstates of any field quadrature operator restricted to a finite 
number TV of photons are studied, in terms of the Hermite polynomials. By (naturally) 
defining approximate eigenstates, which represent highly localized wavefunctions with up 
to TV photons, one can arrive at an appropriate notion of limit for the spectrum of the 
quadrature as TV goes to infinity, in the sense that the limit coincides with the spectrum of 
the infinite-dimensional quadrature operator. In particular, this notion allows the spectra 
of truncated phase operators to tend to the complete unit circle, as one would expect. A 
regular structure for the zeros of the Christoffel-Darboux kernel is also shown. 

1 Introduction 

The optical phase within quantum optics has posed a problem ever since Dirac [1] first for- 
mulated an approach to the subject based on Hamilton's equations for the energy and phase 
as conjugate variables, and using a hermitian operator for the phase which was later proved 
inconsistent. Several approaches to this problem have included the use of unphysical negative 
number states [2], nonunitary phase operators [3,4], and the truncation of the Hilbert space to 
obtain unitary operators. 

This last scheme, due to Pegg and Barnctt [5-7], consists of restricting the full Hilbert space 
H to a subspacc containing at most TV photons, H N = span{|0) , . . . , \N)}, through the action 
of the projector IIn = X^=o l n ) ( n \- Thus a polar decomposition is found for the restricted 
annihilation operator a at = n Aran at, and the unitary exponential phase operator is used to 
make physical predictions (in particular those involving its spectrum and eigenvectors) and 
finally the limit TV — > oo is taken. 

This approach is valuable in that it permits, as opposed to other alternatives, the con- 
struction of a proper Hermitian phase operator, which, however, can be seen as problematic in 
that it depends on a limiting procedure that some authors have found unacceptable [8]. After 
the publication of Pegg and Barnett's work, a number of articles have appeared which treat a 
number of problems related to the truncation of the dimension of the original Hilbert space, 
searching, in particular, for coherent states [9, and references therein]. The relationship of the 
full space with the truncated space, and the commutation relations possible in it, has also been 
studied in [10-17]. 

In this contribution we describe, within this formalism, the spectrum and eigenvectors of 
the position quadrature operator £at for a finite number of photons (and thus for any arbitrary 
quadrature), making extensive use of Hermite polynomials [18, 19]. (These polynomials' con- 
nections with the harmonic oscillator system are deep, and have been studied e.g. in [20,21].) 
The normalization of the eigenstates is resolved using the Christoffel-Darboux kernel, and two 
possible normalizations are given (Section 2). 
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We then show that approximate cigcnstatcs |A)^ can be naturally defined, similarly to 
[18], for all real (pseudo)-eigenvalues A (Section 3), and we define and study a quantitative 
measure, djv(A) = ||(£jv — A) lA)^ || 2 , for the exactness of this approximation. We give suitable 
approximations for this function and find that its analytical continuation into the complex 
plane has poles on the zeros of the Christoffcl-Darboux kernel; we show numerical evidence for 
a regular structure of the latter. 

By studying the zeros of g?at(A), which give the exact eigenstates, one can define an appro- 
priate notion of limit for the spectrum of £tv , in the sense that it tends to the full spectrum of 
the infinite-dimensional position operator £ as A" — > oo. We also show that with this notion of 
limit the finite spectra of the truncated phase operators of Pegg and Barnett (which cover N + 1 
equally spaced points on the unit circle) tend to the complete unit circle - which is what one 
would expect from a phase operator -, in contrast to their assertion in an early paper [22] that 
in the N — > oo limit the phase operator should have only a countable number of eigenvalues. 

A measure of how good an approximation a given vector is to an eigenvector of the position 
operator is also given by the spread of its wavefunction. This is studied in Section 4, where we 
show that the wavefunctions are highly localized at the pseudo-eigenvalues. We also calculate 
and approximate the expectation value and dispersion of £at and £. We finish in Section 5 with 
a few remarks on the Jacobi-matrix representation of the Hermite polynomials. 

This work is presented partly because of its own intrinsic interest, and partly in the hope 
that it will help to understand better the limiting procedure in Pegg and Barnett 's formalism. 
There one tries to understand a hard-to-define operator by approximating it with a sequence 
of finite-dimensional operators; here we have studied the analogous limiting procedure for the 
better-understood position operator. More generally, our results give some indication of features 
to expect when restricting physical operators with continuous spectra to finite-dimensional 
sub-spaces. 

Additionally, there are pedagogical applications to the development of finite-dimensional 
quantum mechanics through the discretization of a particle's position. In that context, the 
constructions in this paper form a more natural basis than that obtained by taking equally- 
spaced points. 

The results presented here constitute an extension to the work in [23] . 



2 Quadrature spectra 



Since an arbitrary quadrature £p = -4= ) is unitarily similar to the position 

quadrature, through the operator U(j3) = e 4,3n , it suffices to consider only the position quadra- 



ture. Thus, we consider the restricted quadrature £jv = ^ 
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This operator is the truncation, through £jy = n^n^y, of the infinite-dimensional position 
operator £, whose domain is £>(£) = E H : (,\ip) £ H}, and whose spectrum is continuous 
and consists of all real numbers, with an infinite-norm eigenbasis |£) normalized to = 

5(Z -O such that £ |0= £10- 
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The characteristic polynomial for £jy (modulo a sign), dct(£ + A), can be expanded by minors 
along the last row and then along the last column to give a recurrence relation, 



2A - 2N 
det(6v+i + A) = — det(&v + A) — -j- det(&v-x + A). 



(2) 



This, together with the initial polynomials dct(£o + A) = A and det(£i + A) = h (2A 2 — l), for 
N = and N = 1 respectively, is sufficient to determine the polynomials for all N. 

This can be used to relate them to the standard families of orthogonal polynomials. Since 
the Hermitc polynomials (as defined in [24]) follow the recurrence relation 



with initial values H\(x) = 2x and H2{x) 
polynomials are given by 

det(&v + A) 



H n+1 {x) = 2xH n (x) - 2nH n _ 1 (x), (3) 
4a; 2 — 2, one can easily prove that the characteristic 
f 



2 N+1 



H 



N+l 



(A). 



(4) 



In particular, this means that the eigenvalues of £jv are exactly the roots of the (N + 1) 
Hermite polynomial. Since there is no exact formula for these roots above the 9th order (until 
which the parity of the polynomials allows exact solutions in terms of radicals) , one might think 
at first that the diagonalization of £jv is an unassailable problem. 

However, the solution can indeed be found by leaving the eigenvalue A alone and reducing 

the system — Xj \X) N = to upper triangular form using Gaussian elimination; after this, 
the system is bidiagonal with coefficients in terms of the Hermitc polynomials. (Here we have 
switched from the combination £jv + A to £jv — A to simplify the algebra, but no results are 
affected since as the Hermite polynomials have definite parity the zeros are unchanged.) The 
solution is then found to be 



y> H n (X) 

2-^ . /on„ ! 



(5) 



in which c/v(A) is a (nonzero) normalization constant which will be discussed later. 

Independently of how it was obtained, expression (5) can easily be seen to be a solution to 

the system f|jv — Xj \X) N = by expressing the latter as the equations 



1|A) 



-A(0|A) w - 
N -X{n\X) N 



V=2 mN 



0. 



(JV - 1|A> „ - A <JV|A) 



N 



(n+l|A) w =0, < n < N, 



= 0. 



(6a) 
(6b) 
(6c) 



The solution (5) satisfies (6a) trivially, and (6b) can be seen to be equivalent to the Hermite 
polynomials' three-term recurrence relation (3). It is only (6c) which depends on A being an 
eigenvalue for £, since from the Hermite recurrence relation 

| (N 1\X) N X (N\X) N = -I^LhwW. 
Thus, for any X the vector |A) ^ differs only in a single component from being an eigenvector 

of i N - 

(&-*)w*~5$& 



H N+1 (X)\N). 



(7) 
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Since the (N + 1) Hcrmitc polynomial has N + 1 zeros, this construction gives as many 
eigenvectors as the dimension of H N and therefore solves the diagonalization problem. The 
normalization, however, is still to be fixed. It is natural to ask that the eigenvectors have unit 
norm: that is, to require that 



CAT 



wi 2 E 



2 n n\ 



n(M^)n 



(8) 



Here the summation can be resolved into a much simpler form through the use- of the Christoffel- 
Darboux formula [24, eq. (22.12.1)] for the Hermite polynomials, 
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E 

n=0 



2 n nl 



1 H N+1 (X)H N (fi) - H N (\)H N+1 (ji) 



2 N + 1 N\ 



A — /i 



which in the limit fx — > A gives 



N 

E 

n=0 



H n {\y _ H N {\)H' N+1 {\) - H N+1 {\)H' N {\) 
2 n nl 2 N + 1 N\ 

= (N + l)gjy(A) 2 - jVg j v +1 (A)g jV _ 1 (A) 
2 N N\ 



(9) 



(10a) 



(10b) 



The Christoffcl-Darboux kernel is of particular importance in this setting, since it is in fact 



part of the integral kernel of the projection operator 11^ = Yln=Q l n ) This can be seen by 

IIjv C ) , using the fact that 



calculating its matrix elements in the position representation, 
the number-basis wavefunctions are given by [25, eq. By. 33] 



^J2 n n\^ 



One obtains 



n 



e) = E «i»> wo = E ^Jfh -^w 



1 H N+1 (QH N (C)-H N (QH N+1 (C) 

2 N + 1 N\y^ ' 

which can be used to obtain the projector's action on wavefunctions as 



-(£ 2 +«' 2 )/2 



n 



N 



2 N + l N\^ 



(11) 

In this context, the normalization constant cat (A) is related to the diagonal matrix elements 
of the projector via 

1 



On 



n 



N 



where A has been replaced by £ in cjy(A). Additionally, if (8) is used as a definition for cat(A), 
one can use this to obtain the identity 



(1A 



MA)| 2 



Tr(njv) = JV + 1. 
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Of course, the normalization jv"(A|A),y = 1 is not the only interesting choice. One other 
important option is to have cjv(A) be independent of TV; in that case, it is particularly fruitful 
to have c^(X) = 7r~ 1 / 4 e~ A / 2 , so that the eigenvectors become 

which arc the truncation through II at of the full (infinite-dimensional) position eigenvectors with 
eigenvalue A. These truncated eigenvectors have a finite norm which is a monotone function 
of TV; in fact, with this normalization, ^(AIA)^ — * 00 as N — > 00, which reflects the fact that 
infinite-dimensional position eigenstates have infinite norm. 



3 Approximate eigenvectors 



For each of the TV + 1 zeros A of ffjv+i.(A) = 0, we have shown the existence of an eigenvector 
\X) N of with £n |A)jy = AIA)^. However, the vectors \X} N can be defined through (5) 
even if A is not an eigenvalue of £/v. In that case \X) N is of course no longer an eigenvector, 
but we have seen that it is "almost" one, in the sense that only one component disrupts this 
behaviour; further, if one is to take the limit TV — > 00 then the |TV) component would seem to 
lose importance. 

Thus, we propose to view the vector |A)jy as an approximate eigenvector for all A, and to 
judge the exactness of the approximation by the function 



d N (\) := |(&r-A) W 



MA)| 2 
2 N + 2 Nl 



H 



N+l 



(A) 2 



(12) 



It is important to note that these are neither eigenvectors nor eigenvalues. The index A is 
perhaps best referred to as a pseudoeigenvalue; while there is a continuum of these, there 
is of course only a finite number of eigenvalues. Analogues of these states have been briefly 
considered for angle variables in [26]. 

It is natural to demand the normalization ^(AIA)^ = 1 so that the length of |A)jy does not 
affect the limiting behaviour of djv(A); one then has 



H 



N+l 



(A) 5 



2 H N (X)H' N+1 (X) - H N+1 (X)H' N (X) ' 
This can be simplified using the identities 



and 



Hn(X)H' n+1 (X) — Hn + i(X)H' n (X) 
H N (X) 1 H' N+1 (X) 



-H N+1 (X) - 



H N (X) 
(A) 



H 



N+l 



1 



Hn+i(X) 
to the remarkable form 



2(TV + l)ff w +i(A) 2(TV + l)dA 
-(TV + 1) 



In Iff 



N+l 



(A)| 



One can also express this as 



dA 2 



^[lnlff^+^A)!]' 



In 



VPWA)i 



(13) 



(14) 
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which immediately gives an asymptotic expression for large A, since then Hj^ + i(X) « 2 N+1 X N+1 . 
After differentiation, this means that c?tv(A) ~ A 2 for large A, which makes sense: one expects 
the vector lA)^ to be an increasingly poor approximation for an eigenvector away from the 
region with eigenvalues. 

However, one can say much more about d/v(A): for instance, expression (13) does not readily 
give much information, and it is not very practical for algebraic manipulation or numerical 
evaluation. For such purposes, it is more convenient to use the expression 



djv(A) 



H 



N+l 



(A) 5 



4 (N + 1)H N (X)H N (X) - NH N+1 (X)H N ^(X) ' 



(15) 



The behaviour of d/v(A) splits into two regions: it is oscillatory for small |A| and quadratically 
growing for large |A|. Specifically, one can prove that 



d N (\) = A — 



37V N(5 - N) 



4A 2 



+ 



d N (X) 



for A > V2N + T, 

l 1 + (-l) N + 1 COS (2V2NT3X (l - 

2 



1 



cos 



V2W+1 



Al 



1 A 2 
12 2W+1 



+ 0{N~ 1 ) 



(16a) 



(16b) 



for A < V2N + T. 

The function djv(A), along with both of these approximations, is shown in figure 1. 

d N (X) 

5 r 

4 : 

3 : 

2 : 
1 





-1 

Figure 1: The function d/v(A) and the approximations (16a), up to the constant term (dashed, 
red), and (16b), with the series truncated as they are shown in the text (dotted, black), for 
N = 15. 

The first of these approximations can easily be expected since d/v(A) is a rational function 
of A; in fact, it is simply a Laurent series, valid for sufficiently large A. As such, it can be 
expressed as djv(A) = J2kLo a k^ 2 ~ 2k , where all the coefficients can be found either manually 
(expressing Hn+i(X) as a sum of A 2 ((N + 1)Hn(X)Hn(\) — NHn+i(X)Hn-i(X)) and terms 
of 0(X 2N ), and so on) or analytically, integrating around a large circle in the complex plane, as 



1 

2ttI 



X 2k ~ z d N {X)AX. 



This reduces in turn to a sum of residues over the poles of c?at(A), which are exactly the zeros of 
the Christoffcl-Darboux kernel, expression (10). Since the kernel is a sum of squares, it is never 
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Figure 2: Zeros of the Christoffcl-Darboux kernel for the Hermite polynomials, given by 
E„=o H n (X) 2 /2 n nl, on the complex A plane, for N = 5, 10, 25, 50, 150 and 500. Care must be 
taken since the scales on both axes do not coincide. 



zero for real variables (a fact which is often used to prove interpolation theorems for zeros of 
orthogonal polynomials) so the poles have nonzero imaginary parts. 

However, despite the extensive attention to the Christoffcl-Darboux kernel in the orthogonal- 
polynomials literature [27], to the best of our knowledge there arc no suitable expressions or 
approximations for these zeros. It is clear, on the other hand, that these zeros do have some 
additional structure, which can easily be seen by plotting them: they appear to lie, with 
increasingly even spacing, on two complex-conjugate curves increasingly close to the real axis 
as N increases. A plot of these zeros for TV = 5, 10, 25, 50, 150 and 500 is shown in figure 2. 

The second approximation, (16b), correct to order 0(N~ 1 ) uniformly for bounded A, is 
obtained from the approximation 

e- x2/2 H n (x) = 2t+3Vd(7T7i)- 1/4 (sin(0))- 1/2 

(sin(2«» - 20) + ^ +0(n- 1 ) 

where x = \J2n + 1 cos(^) and <j> £ [e,7r — e], found in [28, p. 201], by further approximating 
the argument of the sine function using the Taylor series 

4 4 x ( x^ 

sin(20) - 2(j) = -7T + x - — tttttt + O 







/ n 
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V2n+ 1 (27i+l) 3 /2 6 y n 5/2 



to obtain 



2/2 H n (x) 



2 n n!v / 2 / 2n+ 1 



1/7TO \ 2n + 1 — x 



x cos V271 + lx 1 - 



1/4 



1 X 2 A 7T 

n 



6 2n + 1 J 2 



One consequence is that the spacing of the zeros of djsr(X) (which give the eigenvalues of £n) 
decreases as tt/ \/2N + 1 for large N, and that this spacing is increasingly even as N increases. 
Furthermore, the region with zeros covers the increasingly large intervals —\/2N + 1 < A < 

vWTT. 

One can then show that, given any real number A, a positive number e and an integer No, 
there exists an eigenvalue Ao of £at with N > No such that |A — Ao| < e. This can in turn be 
used to define an appropriate notion of the limit of the spectra c(£jv), with the desired property 
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that the spectra tend to the spectrum of the full, infinite-dimensional, position operator £: 



lim a(i N ) := f) U ^ = R = a(£), (17) 

fc=0 N=k 

where the overline means topological closure in K. 

This notion of a limit of a sequence of sets means that the points in the limit set are 
arbitrarily close to sets of the sequence with arbitrarily large indexes, and it resembles the limit 
superior [29, p. 16] of a sequence of sets, which omits the topological closure. Additional uses 
for this definition can be easily found: for example it is in this sense that the hypcrboloids 
x 2 + y 2 — z 2 = ±a 2 "tend to" the cone x 2 + y 2 = z 2 as the waist a tends to zero, and that the 
ellipsoids x 2 + y 2 + z 2 /c 2 = 1, when c — > oo and they become increasingly elongated, "tend to" 
the cylinder x 2 + y 2 = 1. 

Our definition (17) also allows the finite spectra of the truncated phase operators of Pcgg and 
Barnctt (which cover N + 1 equally spaced points on the unit circle) to tend to the complete 
unit circle, which is what one would expect from a phase operator. This is in contrast to 
Pegg and Barnett's assertion [22] that in the N — s- oo limit the phase operator should have 
only a countable number of eigenvalues: it is clearly dependent on the notion of limit used, 
which should reproduce in the simpler example of the truncated quadratures £at the required 
behaviour. 

Another limit that deserves notice is that of \X) N as |A| 3> y/2N + 1, where the number-basis 
components in (5) corresponding to \n) with n < N, which scale as 2" A™, are drowned out by 
the \N) contribution. One therefore has, up to a phase, lA)^ — > \N) as A — > ±oo. This limit 
will become apparent on the graphs of the wavefunctions (ClA)^ in figure 3. 

One further detail of the asymptotic behaviour of d^(X) is noteworthy, and it relates to the 
amplitude of the oscillations, which is basically independent of N: div(A) oscillates sinusoidally 
from to 1, instead of - as one might hope - tending uniformly to zero. This means that for 
any real A, there will always be values Ai arbitrarily close to A that will have cIn (Ai) = 1, with 
arbitrarily large TV. 

The fact that d/v(Ai) = 1 = ^(AilAi)^ can be seen as implying that |Ai)jy is as far away 
as it can be from being an eigenvector, which would disqualify the curve An- | A) ^ in H N from 
describing approximate eigenvectors. However, this overlooks the fact that £n has a relatively 
large (operator) norm, and that one can find states such as \N), the last number-basis vector, 

for which 



£ — A) |7V) = A 2 + y. Since \N) represents the wavefunction within % N with 



2 

2 

arguably the largest position dispersion, it can be seen that the vectors lA)^ are in fact rather 
close to being position eigenvectors. 



4 Wavefunctions and moments 

As we mentioned at the end of last section, one qualitative measure of how good an approxima- 
tion to an eigenvector of £jv is a given vector is the spread of its wavefunction. In this respect, 
one might hope that the vectors lA)^ represent wavefunctions highly localized around A with 
small dispersion. 

The wavefunctions can be found, again using the Christoffcl-Darboux sum formula, now in 
its bivariate form, (9): 

,, m m f ^(A),,, M ^ g n (A)g»(Q e-?l 2 
{i\\) N - c N (X) ^ -= (£|n) = c N (X) ^ — 

n— ' n— 

Cjv (A) e-?l 2 H N+1 (\)H N (£) - H N (X)H N+1 (£) 

2 N+i m n l/4 X -£ 
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If one chooses cjy(A) to be real and positive then one has 

e-« 2/2 1 H N+1 (X)H N (0 -H N (X)H N+1 (0 



^W+mV^ A - £ ^ Hn{X )H' n+1 (\) - H' N (X)H N+1 (X)) ' 



(18) 



The graphs of selected wavefunctions are shown in figure 3; however, it is more instructing to 
see the wavefunctions change with A in a continuous animation (see supporting material online). 



<f|A) A 



<f|A) 



1(1 5 -10 





(b) A = 
<f|A>« 




(d) A = 2 




(c) A = 5 (f) A = 10 

Figure 3: Wavefunctions (£|A)at as functions of £ for N = 16 and the indicated values of A. 

Seeing now that the wavefunctions are highly localized, it is interesting to calculate the 
expected value of the position and its dispersion. There is an initial question of whether to use 
the truncated position, £jv, or its full version, £, but it turns out that their expected values 
coincide: since \X) N is in H N , the range of II tv, 



N 



X) = (X 

N N 



n^n 



A) = (A 

N N 



N 



N 



This is convenient, since we know already from (7) the action of £at on the vectors \X} N . We 
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thus have 



N 



IN 



N 



N 



= X -2 



1 cat (A) 
l|c w (A)| 2 



1 c N (X) 



H 



N+l 



2 N Nl 



H N+1 (X) ■ N (X\N) 

Hn+i(X)Hn(A) 
Hn+i{A)Hn(\) 



(A) \N) 



2 {N + l)H N (\) 2 - NH N+ i(\)H N -i(\) ' 

In particular, the expected value (£jv) coincides with the pseudo-eigenvalue A not only on 
the zeros of iZjv+i! when it should since |A)jy is then an eigenvector, but also on the zeros of 
fljv- This is important since these second values are interspersed with the first and therefore 
occur on the peaks of the oscillations of djv(A), a further indication that the treatment of \X) N 

as an approximate eigenvector is appropriate. 

Approximating the expression for the expected value of the position using the same methods 
that were used to approximate (ijv(A) by expressions (16), one obtains 



N(N - 2) 



+ ■ 



V 2A 3 
for A > V2N + 1, and 

A 



(0 « A - 



(A( 1 + safe + -")) + (- 1 )' fain ( 2 ^ FFlA ( 1 + ^ + '")) 



^VTT [l + cos ( 7 =| TT A (l + i g^ij + • ■ • ))] 



for A < v / 2JV+T. 



(19a) 



(19b) 



The graph of (A 



A' 



A ) as a function of A is shown, along with both of the above aproxi- 

N 



mations (with the series truncated as shown above), in figure 4. 




Figure 4: Expected value (£jv) as a function of A for = 16, with the approximations (19a) 
and (19b) shown dashed and dotted, respectively 

Since the amplitude of the oscillations in (19b) does decrease uniformly, in the A <§C y/2N + 1 
regime one can consider A to be an excellent approximation to the expectation value (£jv)- In 
general, jv(A|£jv|A)iV covers all of the values between the smallest and largest eigenvalue of £jv, 
i.e. the first and the last zeros of i/jv+i (which is approximately the region —\J2N + 1 < A < 
\/2N + 1) linearly in A except for oscillations whose amplitude tends to zero, although near the 
edges of this interval it can occur that more than one A correspond to the same expectation 
value, which happens because the amplitude of the oscillations increases with A. 
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The decay of (£jv) to zero as A — > oo is also interesting in its own right, since it represents 
the fact that in that regime the states default to the last number ket, \N), which has expected 
value for the position. 

We close this section by giving an explicit expression for the position dispersion. This is 
complicated by the fact that the second moments of £ and £jv do not coincide, since Il7vC 2 nw 7^ 
In fact, one has 



A) = (A 

N N 



n w £n^n w 



A + A 



/v 



A? 



N 



N 



C2 



N + 1 9 



' N 



We therefore calculate 



N 



X) = (X 

N N 



A 



2 1 \cnW\* 



2 2 N Nl 

which gives, for the position uncertainty, 



H N+1 {X) {XH N (X) + NHn-^X)) , 



N 













(liv-(|iv>) 2 


a\ = (A 


a 


A) - (a 




In nX 


IN N\ 





N 



MA) ( 1 - %i!^(A) 2 



2 N Nl 



This can be approximated by the expressions 




Figure 5: Position dispersion (A£at) 2 as a function of A for N = 8, along with the (dashed and 
dotted) approximations (20). 




N N(N + 3) 
4A 1 4 

d N (X)+0(N~ 1 ) 



for A > V2JV+1, and (20a) 
for A < V2N + 1, (20b) 



11 



On the spectrum of field quadratures for a finite number of photons 



which are shown in figure 5. 

The full position dispersion, (j^^j , is slightly more complicated. As said above, it differs 
in one term from the truncated-position dispersion, and can therefore be written 



iV 



£-<0 



d N (X) 1 



iY 



MA)P 

2 N N\ 



H N (\)' 



N+l\c N (\)\' 
2 2 N N\ 



(21) 



Because of the N + 1 factor, the last term does not tend to zero as N — > oo and is therefore 
not negligible. Furthermore, since the extra term is proportional to Hn(\) and the original 
cf/v(A) is proportional to iJjv+i(A), in the oscillatory region the two contributions are out of 
phase and therefore interfere constructively to make an (almost) constant function close to 
A = 0. As A grows, the differing frequencies of oscillation of both terms make them interfere 

and (j^^j docs show some oscillations. This can be seen in figure 6, and is made precise by 
the approximations 

2N + 1 N(N + 2) N(2N 2 + 27V - 9) 



(MY 



(AO 2 



2A 2 



4A 4 



for A > V2N + T, and 
1 + {~1) N sin (2y/2N + 2X (l 



6 2JV+2 



•))' 



V2AT+2 



1 + 



6 2JV+2 



0) 



. V2AT+1 



All 



12 2AT+1 



■)) 



for A < V2N + 1. 



(22a) 



(22b) 



(A£) 2 




Figure 6: Dispersion of the (non-truncated) position, (A£) 2 , for TV = 16, together with the 
approximations (22a), dotted, and (22b), dashed. At A = the dispersion is |, independently 
of N, and there the function is increasingly flatter as N — > oo. 



For large A, both dispersions are bounded, which they must since they are the expectation 
values of finite-dimensional operators. It is interesting to note that the full-position operator 
dispersion, (A£) 2 , seems to be bounded from below, in the sense that (A£) 2 — s- | uniformly in 
any compact set A G [— Ao,Aq], which contrasts with the fact that £ has eigenstates with zero 
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dispersion. Although the latter do have infinite norm, it is well known that squeezed states 
with finite norm and position dispersion below the Hciscnberg limit of ^ exist. 

For completeness, we include here some additional results. The inner product between 
different (pscudo-)eigenvectors \X) N is given by 



Cj v(A')* CjV (A) H N+1 (\)H N (\>) - H N (\)H N+1 (\>) 
M A Wn - 2 N + 1 N\ ~\~\' ' 



in particular, it is zero when both A and A' are roots of Hn + \. but also when they are roots 
of Hm, which is a consequence of the fact that they are eigenstates of £n-i, and are therefore 
well treated as pseudo-eigenstates of £jv- The matrix element 



'A' 

N 



_ c N {\')*c N {X) \'H N+1 {\)H N {\') - \H N {\)H N+l {\') 
>n~ A -A' [ ' 



is also of interest; the matrix element ^(A'lTrlA)^ can also be calculated 



5 The Jacobi matrix 



In this section we briefly explore some mathematical results related to the work described in 
this article. Most importantly, we retake the expression for the characteristic polynomials of 
the truncated position, £jv, given in (4), and which can be written as 



/ A 4= 



H N+1 (X) = 2 N+1 dctfev + A) = 2 N+1 det 



V 



1 

^2 





1 
A 





V2 

A 



\ 





^ A 

\/2 



(25) 



This can now be turned around, and viewed as a representation of the Hcrmitc polynomials. 
This representation in terms of the determinant of a tridiagonal matrix exists for all families 
of orthogonal polynomials and is known in the special-functions literature as the Jacobi-matrix 
representation [30]. Jacobi himself used this expression, as well as the explicit diagonalization 
given above of the tridiagonal matrix, to study Gaussian quadrature [31], that is, the approxi- 
mation of the integral of a function by a finite sum of evaluation on a number of points. These 
points are the zeros of appropriately chosen polynomials, which can be efficiently evaluated by 
this method since it uses the algorithms of linear algebra. 

An interesting application of this development comes if one couples the expression (25) with 
the Caylcy-Hamilton theorem, which states that every matrix satisfies its own characteristic 
polynomial. This means that the truncated quadratures £jv arc matrix solutions to the equation 



H N+1 (X) = 0, 



(26) 



which is not trivial, particularly when one considers that for large N the real solutions are 
known to exist but there is no known simple or elementary exact formula for them. 

Furthermore, since each zero of Hn+i appears exactly once as an eigenvalue of ^n, this 
means that the Hermite polynomials are not only the characteristic polynomials of the matrices 
£at, but they are also their minimal polynomials (the minimal polynomial p of a matrix A is 
the unique monic polynomial of minimal degree for which p(A) = 0). Although there are other 
matrix zeros of (26), the condition of having the Hermite polynomials as minimal polynomial 
specifies the truncated quadratures £at uniquely, up to matrix equivalence. 
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Finally, the work in this article may be used to obtain results for the Hcrmite polynomials 
which would be hard to obtain in other ways. We give one example: by numbering the zeros A& 
of Hn + i(\) = with k = 0, 1, . . . , N we obtain an othornormal basis {|Afc) N : k = 0, 1, . . . , N}, 
in terms of which the orthonormalization of the number kcts reads 



N < N 



S mn = (m\n) = V (m|A fe )jv N (Xk\n) = -==== V \c N (X k )\ 2 H m (X k )H n {X k ), 
, n \J2 m ra\2 n n\ , n 

k—0 k—0 

an expression in which the sum is over all the zeros of Hm+i and would therefore be particularly 
hard to obtain by ordinary means. 



6 Conclusions 

We have obtained the spectrum of the position quadrature operator £jy for a finite number 
N of photons in terms of the zeroes of Hermite polynomials, as well as its eigenvectors in 
terms of lower-degree Hcrmite polynomials evaluated on these zeros. The normalization of the 
cigenstates is given through the Christoffcl-Darboux kernel, and a regular structure for these 
is shown numerically. Approximate eigenstates for £jy were naturally defined, which represent 
highly localized wavefunctions centered around any value between the least and greatest zeros 
of Hn + i. By using an appropriate notion of limit for these spectra, in the sense that they 
tend to the full spectrum of the infinite-dimensional position operator £, we showed that the 
finite spectrum of the truncated phase operators tends to the complete unit circle, as one would 
expect. The same procedure can be followed for any arbitrary field quadrature. 



Acknowledgments 

The authors wish to thank A. Frank and K. B. Wolf for discussions. This work was partially 
supported by CONACYT, SEP and by DGAPA-UNAM under project IN102811. 



References 

[1] P A M Dirac. The quantum theory of the emission and absorption of radiation. Proc. R. 
Soc. Lond. A 114 (1927), pp. 243-65. 

[2] S M Barnett and D T Pegg. Phase in quantum optics. J. Phys. A 19 (1986), pp. 3849-62. 
doi: 10. 1088 /0305-4470/19/ 18/030. 

[3] L Susskind and J Glogower. Quantum mechanical phase and time operator. Physics 1 
(1964), pp. 49-64. 

[4] P Carruthers and M M Nieto. Phase and angle variables in quantum mechanics. Rev. 
Mod. Phys. 40 (1968), pp. 411-440. doi:10.1103/RevModPhys.40.411. 

[5] M M Nieto. Quantum phase and quantum phase operators: some physics and some 
history. Phys. Scripta T48 (1993), pp. 5-12. doi:10.1088/0031-8949/1993/T48/001, 
arXiv:hcp-th/9304036vl. 

[6] D T Pegg and S M Barnett. Tutorial review: Quantum optical phase. J. Mod. Opt. 44 
(1997), pp. 225-64. doi:10. 1080/09500349708241868. 

[7] S M Barnett and J A Vaccaro. The quantum phase operator: a review. Series in optics 
and optoelectronics. Taylor & Francis, New York, 2007. 



14 



E PlSANTY AND E NAHMAD-ACHAR 

[8] M A Hcnnings, T B Smith, and D A Dubin. Approximations to the quantum phase opera- 
tor. J. Phys. A: Math. Gen. 28 (1995), pp. 6809-6856. doi:10.1088/0305-4470/28/23/026. 

[9] B Roy and P Roy. Coherent states, even and odd coherent states in a finite-dimensional 
Hilbert space and their properties. J. Phys. A: Math. Gen. 31 (1998), pp. 1307-17. 
doi:10. 1088/0305-4470/31/4/018. 

[10] S M Barnett and D T Pegg. Limiting procedures for the optical phase operator. J. Mod. 
Opt. 39 (1992), pp. 2121-29. doi:10.1080/09500349214552141. 

[11] D T Pegg, J A Vaccaro, and S M Barnett. Quantum-optical phase and canonical conjuga- 
tion. J. Mod. Opt. 37 (1990), pp. 1703-10. doi:10.1080/09500349014551931. 

[12] W C Hammel. Finite canonical commutation relations, 1999. Available online at 
http : / / graham . main . nc . us/~bhammel/FCCR/f ccr . html#contents. 

[13] T Santhanam and A Tckumalla. Quantum mechanics in finite dimensions. Found. Phys. 
6 (1976), pp. 583-587, 1976. doi:10.1007/BF00715110. 

[14] S Gudder and V Naroditsky. Finite-dimensional quantum mechanics. Int. J. Theor. Phys. 
20 (1981), pp. 619-643. doi:10.1007/BF00671374. 

[15] A C de la Torre and D Goyeneche. Quantum mechanics in finite-dimensional hilbert 
space. Am. J. Phys. 71 (2003), pp. 49-54, 2003. doi:10. 1119/1. 1514208, arXiv:quant-ph/ 
0205159vl. 

[16] N Cotfas, J P Gazeau, and A Vourdas. Finite-dimensional hilbert space and frame quantiza- 
tion. J. Phys. A: Math. Theo. 44 (2011), p. 175303. doi:10.1088/1751-8113/44/17/175303. 

[17] R Jagannathan, T S Santhanam, and R Vasudevan. Finite-dimensional quantum mechanics 
of a particle. Int. J. Theor. Phys. 20 (1981), pp. 755-773. doi:10.1007/BF00674253. 

[18] H A Buchdahl. Concerning a kind of truncated quantized linear harmonic oscillator. Am. 
J. Phys. 35 (1967), pp. 210-18. doi:10. 1119/1. 1974004. 

[19] P Figurny, A Orlowski, and K Wodkiewicz. Squeezed fluctuations of truncated photon 
operators. Phys. Rev. A 47 (1993), pp. 5151-57. doi:10.1103/PhysRcvA.47.5151. 

[20] V V Dodonov, O V Man'ko, and V I Man'ko. Photon distribution for one-mode mixed 
light with a generic Gaussian Wigner function. Phys. Rev. A 49 (1994), pp. 2993-3001. 
doi:10.1103/PhysRevA.49.2993, arXiv:hcp-th/9308004v2. 

[21] V V Dodonov, O V Man'ko, and V I Man'ko. Multidimensional Hcrmite polynomials 
and photon distribution for polymode mixed light. Phys. Rev. A 50 (1994), pp. 813-17. 
doi:10.1103/PhysRevA.50.813, arXiv:hep-th/9308033v2. 

[22] D T Pegg and S M Barnett. Unitary phase operator in quantum mechanics. Europhys. 
Lett. 6 (1988), pp. 483-87. doi:10. 1209/0295-5075/6/6/002. 

[23] E Pisanty. Estados coherentes generalizados y estructura analftica del operador de aniqui- 
lacion. Bachelor's thesis, Universidad Nacional Autonoma de Mexico, Mexico, 2011. Avail- 
able online from http : //132 . 248 . 9 . 195/ptb2011/septiembre/0673328/Index .html. 

[24] M Abramowitz and I A Stegun, editors. Handbook of Mathematical Functions, with For- 
mulas, Graphs, and Mathematical Tables, ninth printing. Dover, 1972. 

[25] C Cohen- Tannoudji, B Diu, and F Laloe. Quantum mechanics. Wiley, New York, 2 nd 
edition, 1977. 



15 



On the spectrum of field quadratures for a finite number of photons 



[26] S M Barnctt and D T Pcgg. Quantum theory of rotation angles. Phys. Rev. A, 41 (1990), 
pp. 3427-3435. doi: 10.1 103/PhysRevA.41. 3427. 

[27] B Simon. The Christoffel-Darboux kernel, in Perspectives in partial differential equations, 
harmonic analysis and applications, volume 79 of Proc. Sympos. Pure Math., pp. 295-335. 
American Mathematical Society, Providence, RI, 2008. 

[28] G Szego. Orthogonal Polynomials, volume XXIII of Colloquium Publications. American 
Mathematical Society, 4 th edition, 1975. 

[29] P Halmos. Measure theory. Springer- Ver lag, New York, 1974. 

[30] G E Andrews, R Askcy, and R Roy. Special Functions, volume 71 of Encyclopedia of 
Mathematics and its Applications. Cambridge University Press, 1999. 

[31] W Gautschi. Orthogonal Polynomials, Computation and Approximation. Numerical Mathe- 
matics and Scientific Computation. Oxford University Press, 1 st edition, 2004. 



10 



